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Abstract 

Using the Green's function associated with the one-dimensional 
Schroedinger equation it is possible to establish a hierarchy of sum 
rules involving the eigenvalues of confining potentials which have only 
a boundstate spectrum. For some potentials the sum rules could lead 
to divergences. It is shown that when this happens it is possible to 
examine the separate sum rules satisfied by the even and odd eigen- 
states of a symmetric confining potential and by subtraction cancel 
the divergences exactly and produce a new sum rule which is free of 
divergences. The procedure is illustrated by considering symmetric 
power law potentials and the use of several examples. One of the ex- 
amples considered shows that the zeros of the Airy function and its 
derivative obey a sum rule and this sum rule is verified. It is also 
shown how the procedure may be generalised to establish sum rules 
for arbitrary symmetric confining potentials. 

PACS: 02.30.Gp, 03.65.Fd, 03.65-w, ll.30.Pb 
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1 Introduction 



Green's functions G are used to study the solutions to inhomogeneous differ- 
ential equations which satisfy homogeneous boundary conditions. G satisfies 
the same boundary conditions as the solution ^ of a second order diffren- 
tial equation. G satisfies a corresponding differential equation with a delta 
function source term (Morse and Fcshbach (1953)). For one-dimensional sec- 
ond order differential equations there are two methods for constructing the 
Green's function. In the first method G is constructed as a sum involving the 
eigenfunctions and eigenvalues of the homogeneous differential equation. In 
the second method G is constructed from two solutions of the homogeneous 
differential equation satisfying appropriate boundary conditions. The exis- 
tence of two alternative ways of constructing the Green's function leads to 
the establishment of certain equalities involving the eigenfunctions and eigen- 
values of the differential operator. For studying the boundstates of quarks 
with definite angular momenta the radial Schroedinger equation with a con- 
fining potential has been used extensively. Such solutions for zero orbital 
angular momentum [viz) the s-states, can be viewed as the odd states of a 
symmetric confining potential in the infinite space — oo < a; < oo (Quigg et 
al (1980)). Since confining potentials have no scattering states and have only 
boundstates the corresponding Green's functions lead to a hierarchy of sum 
rules which only involve the boundstate eigenvalues. Such sum rules can be 
powerful tools to establish the extent to which a set of numerically computed 
eigenvalues approach completeness. A leading order sum rule has been used 
recently by Mezincescu (2000) to shed light on the hypothesis advanced by 
Bender and Boettcher (2001) that a cubic oscillator with imaginary coupling 
strength has a completely real spectrum. 

The sum rules involving all the boundstate eigenvalues do not always 
lead to convergent expressions and prove useful. When this happens it is 
still possible to extract something useful from the Green's functions. In this 
paper it is shown that it is possible to establish separate sum rules for the odd 
and even states in a symmetric potential by considering Green's functions 
satisfying different boundary conditions and that when divergences in the 
sums arise they may be of the same type for the odd and even states. If 
this is the case then it is possible to subtract the two divergent sum rules 
to establish a convergent sum rule involving both the odd and even states. 
The plan of this paper is as follows: the construction of the Green's function 
and the derivation of the sum rules for the inverses of the eigenvalues of 
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confining potentials are presented in Section 2. The sum rules for the even 
and odd states in symmetric confining potentials of the power law form are 
constructed in Section 3. The spectral zeta function for power law potentials 
has been studied by Voros (2000). In Section 4 several solvable examples 
are considered and the sum rules are verified. In Section 5 the case of an 
arbitrary symmetric confining potential is considered and general results are 
presented. 

2 Green's Functions and Sum Rules 

Consider a second order differential operator of the form 

^(^) = -^ + ^(^) (1) 
with a complete set of solutions of the eigenvalue equation 

L^n = A„^„ (2) 

for a choice of two boundary conditions Bi and B2 on ^' at Xi and X2 respec- 
tively. Bi and B2 can be chosen to be homogeneous boundary conditions 
of the Neumann or Dirichlet kind, i.e., the eigenfunction or its derivative 
vanishes at Xi or X2- The orthonormality of the eigenfunctions is given 

by 

/ (x) (x) = 5mn- (3) 

J Xl 

A Green's function G satisfying the same boundary conditions Bi and B2 as 
^„ at Xl and X2 , respectively, may be defined as a solution of 

{L{x)-\)G{x,y;X) = S{x-y). (4) 

A representation of G in terms of the complete set of eigenfunctions and 
eigenvalues is discusssed in standard textbooks on Green's functions. For 
a general V (x) the complete set of eigenfunctions include eigenstates for 
both discrete and continuous eigenvalues. However, if V (x) is a confining 
potential, then the complete set is made up of only discrete eigenvalues A„ 
and G may be represented as a discrete sum given by 

G{x,y;\) = — — " ^5) 
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Using the completeness relation 

Y,'^n{x)^niy) = H^-y) (6) 

n 

and eq.(2) it is easy to verify that G defined by eq.(5) satisfies eq.(4). 

For one-dimensional equations such as eq.(4) an alternative procedure 
to construct G is to solve eq.(4) directly. Let 0i and 02 be two linearly 
independent solutions of 

(L - A) 01,2 = (7) 

chosen such that 0i satisfies the boundary condition Bi but not B2 and 02 
satisfies B2 but not Bi. Then in terms of 0i and 02, G may be constructed 
as 

r^f^...\\ 01 ix<A) 02 (a:>,A) 

'""■^■^' = mK^k)— ''^ 

where the Wronskian W given by 

w^^ ^^ ^ ^ \x#2(a;,A) c^0i(a;,A) 

(01, 02) = 01 (a:, A) 02 (a;, A) — (9) 

satisfies 

dW , , 

and a:< (x>) is the smaller (larger) of (x,y). Using eqs.(7) and (10) it is 
easy to verify that G defined by cq.(8) is a solution of eq.(4). When the 
Green's function can be constructed from eq.(8) then eq.(8) may be taken 
to be the definition of G and eq.(5) may be viewed as the expression of an 
equality relating G to the eigenf unctions and eigenvalues. The existence of 
two expressions for G leads to certain equalities. Using the orthonormality 
condition defined by eq.(3) it can be shown from eq.(5) that 

G{x,x-\)dx = Y,- -. (11) 



Eq.(ll) is an expression of a sum rule if eq.(8) for G is used to evaluate 
the integral. The above procedure can be generalised to provide a hierarchy 
of sum rules as discussed by Sukumar (1990). The orthonormality of the 
eigenfunctions can be shown to lead to a second order sum rule in the form 

/•X2 rX2 -1 

/ dx dyG{x,y;X)G{y,x;X) = J2j: ^72- (12) 

Jx^ Jx^ „ (An — A) 
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Similar sum rules involving a chain of n Green's functions starting from x 

and returning to x via (n — 1) intermediate set of points may be used to 
eatablisli sum rules of arbitrary order n. The sums over the powers of the 
inverses of the eigenvalues for the case A = are referred to as generalised 
zeta functions in the literature. All the sum rules defined so far converge for 
finite values of Xi and X2. For general confining potentials with xi — oo 
and X2 ^ oo extra convergence conditions are needed to ensure that the sum 
rules lead to convergent expressions. These questions are examined in the 
next section. 



3 Sum Rules for Power Law Potentials 

Consider symmetric potentials of the power law form The eigenstates 

split into groups of odd states satisfying the boundary conditions 

Lf^^o* (a;) ^ , Li^^oo* {x) ^0 (13) 
and even states satisfying the boundary conditions 

Lt,_o^^ (x) ^ , Lt^^oo^ (x) ^ . (14) 

If the parameter A is chosen to be equal to then the zero energy Schroedinger 
equation 



dx'^ 

under the substitutions 



* = (15) 



2 1 1 

V = — — - ^ z^vx^ , 13^ — — - , * = (16) 

transforms into 

which is the differential equation satisfied by the modified Bessel functions ( 
Abramowitz and Stegun (1965)). Hence the solutions are of the form 

* =y/x I±f3 [i^x^^ (18) 
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from which we can construct solutions which satisfy appropriate boundary 
conditions as a; — and x oo. The solution which shows decaying be- 
haviour as X — s> oo is \fxKfi i^z) , the solution which vanishes as x ^ is 
-Jxlj^ji (z) and the solution for which the derivative vanishes as x — > is 
y/xl-i3{z). Hence the Green's function appropriate for the description of 
eigenstates satisfying boundary conditions analogous to those in eq. (13) is 

Gi{x,y) = v^^I+p(vx^ K0(pxl^ (19) 

and the Green's function appropriate for the description of eigenstates satis- 
fying boundary conditions analogous to those in eq.(14) is 

G2{x,y) = u^^I_p(px^ Kp(uxi^ . (20) 

The boundary conditions satisfied by Gi and G2 are 

Ltx^o Gi {x, y) ^ , Lt Gi {x,y) 
LU^o -^G2{x,y) ^ 0, LU^^G2{x,y) ^ 0. (21) 

The WKB result for the eigenvalues A„ in the hmit of large quantum number 
n can be derived from the Bohr-Sommerfield formula and is of the form 




(22) 

which shows that the sums over the inverses of the eigenvalues will not con- 
verge if < 2. This also indicates that the integrals of the Green's functions 
over the infinite domain [0, 00] will not converge ii N < 2. When the Green's 
functions in eqs.(19) and (20) are used in eq.(ll) with the choice of limits 
Xi = and X2 = 00 the resulting integrals can be performed (Gradshteyn 
and Ryzhik (1965)) ii N > 2 (i.e) P < ^. The resulting sum rules are 



for the odd states and 



c _ 1 _ f,2-,p r(2/j)r(/j)r (1-4/3) i 

^'"^A^;-^ r (1 - 3/3) r (1 - 2/3) ' ^^^^ 



n=0 
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for the even states. The WKB approximation for the eigenvalues given by 
eq.(22) indicates that when N < 2 the sums given by 5*1 and 5*2 diverge 
but it also indicates that if we look at the difference 5*2 — 5*1 the terms in 

_ 3JV+2 

the resultant sum have the large n behaviour of n ^+2 which would lead to 
convergent sums for all positive definite values of N > 0. Hence we examine 

/"OO 

^2-^1= / {G2{x,x) - Gi{x,x))dx . (25) 
Jo 

Using the addition formulae satisfied by the modified Bessel functions which 
appear in the integrand the integral can be written as 

4/3 ( 1 \ 

S — S2 — S\ = — simrf] I X K"^ [I'x'^j dx (26) 

and the resulting integral can be performed (Gradshteyn and Ryzhik (1965)) 
to obtain the result 

n=oo ^_^„ ^ r (3/3) (2/3)r (/3) sinvr/^ 



n=0 



By using the reflection property satisfied by the gamma functions 

V{z) V{l-z) = (28) 

Sm TTZ 

the results for the sum rules can be simplified. S can be represented in the 
form 

hi r(4/3)r(i-/3) ' f <2- ^""i 

In terms of S the other sum rules derived in this section may be written in 
the compact form 

S2 = y^ =^^^-,iV>2, f/3<iV (30) 
^ ^2n sinnP 2cos27r/3 ' ' V ^' 



Si ^ 1 = S ^ , N>2 , [P<-) . (31) 

^ A2„+i 2cos27r/? V 4; ^ ^ 

We have thus provided a direct proof of the sum rules listed by Voros (2000) 
by a method that is capable of being generalised to other symmetric confining 
potentials. In the next section special values of the power law index for 
which the eigenvalues are readily available in the literature are considered 
and the resulting sum rules are examined. 
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4 Examples of Sum Rules 



4.1 Linear Potential V(x)=x 

For the linear potential the Schroedinger equation for zero energy is given by 
eq.(15) with = 1 and can be identified as the differential equation satisfied 
by Airy functions (Abramowitz and Stegun (1965)). The eigenvalues of the 
even eigenstates whose derivatives vanish at x = are given by the negative 
of the zeros of the derivative of the Airy function and the eigenvalues of 
the odd eigenstates whose wavcfunctions vanish at x = arc given by the 
negative of the zeros of the Airy function. The zeros of the Airy function 
and its derivative are listed in Abramowitz and Stegun (1965) and the first 
10 zeros are given below in Table 1. 



n 


^2n 


A2n+1 


^2n ^2n+l 





1.01879 


2.33811 


0.55386 


1 


3.24820 


4.08795 


0.06324 


2 


4.82010 


5.52056 


0.02632 


3 


6.16331 


6.78671 


0.01490 


4 


7.37218 


7.91113 


0.00977 


5 


8.48849 


9.02265 


0.00697 


6 


9.53545 


10.04017 


0.00527 


7 


10.52766 


11.00852 


0.00415 


8 


11.47506 


11.93602 


0.00337 


9 


12.38479 


12.82878 


0.00279 



Table 1: The eigenvalues of the lowest 10 even and odd eigenstates for the 
linear potential and the differences between the inverses of the eigenvalues 
taken in pairs 

The eigenvalues of the states with larger values of n can be accurately 
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estimated from the asympototic formulae 



2 2 

Lt^^ooX,n = (^) (1 + 4;^) (32) 



2 2 
3 



(Stto \ ^ / 3 

which show that the sums Si and S2 diverge. However if we look at the dif- 
ference S between the sums over the inverses of the even and odd eigenvalues 
the terms in S have the asymptotic behaviour 

it»-«(Aj--A-,) = (J;^)'3^ (34) 

which indicates that the integral over n would converge. The first 10 terms 
in the scries for 5" are tabulated in the last column of Table 1 and add up 
to ~0.691. The remaining terms in the sum can be added by converting the 
sum over n into an integral over n and using the asymptotic estimate given 
by eq. (34) for the integrand. Such an evaluation gives 

f°° 1 / 2 \ ^ 

R = I - ^-0 * = 2 (^j ■ (35) 

Hence an estimate of S is given by choosing n — 10.5 in eq.(35) which gives 
the result that 

00 

S i - ^2n+i ) ~ -691 + 0.037 = 0.728 (36) 

n=Q 

in good agreement with the exact result derivable from eq.(29) using the 
value /3 = I for = 1 which gives 

S = (-Y ^Iv ~ 0.729 . (37) 

\3j r(|) ^ ' 

By calculating more zeros of the Airy function and its derivative exactly 
the sum involving the eigenvalues can be evaluated more accurately without 
approximating part of the sum by an integral over n and the agreement 
between the two sides of eq.(29) can be verified to more significant figures. 
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The asymptotic limits given by eqs.(32) and (33) also indicate that even- 
though the sums over the inverses of the eigenvalues will not converge they 
also show that the sums over squares and higher powers of the inverse of the 
eigenvalues will converge. Hence higher order sum rules such as those arising 
from eq.(12) may be verified by using the appropriate eigenvalues in the sum 
and the appropriate Green's functions in the integral. 

4.2 Simple Harmonic Oscillator V(x) = + 1 

For the oscillator potential V (x) — x^ the eigenvalues are given by A„ = 
2n + 1. The sums 5*1 and 5*2 are both divergent but the difference S gives 

S-.-\^\-\. .... = \ (38, 

in agreement with the value from eq.(29) for /3 = | corresponding to A'^ = 2. 
The zero energy solutions to the Schroedinger equation have to be calculated 
numerically to verify the integral involving the Green's functions. However 
analytical results are possible for the oscillator shifted in energy by one unit. 
Therefore we consider this case in detail. For the shifted simple harmonic 
oscillator potential V [x] = a;^ + 1 the eigenvalues are A„ = 2n + 2 with n 
taking positive even and odd values. The sums S2 and 5*1 over the inverses 
of the even and odd eigenvalues are both logarithmically divergent. However 
the difference between the sums gives 

^ i-f 1111 ln2 

n=0 " 

The zero energy solutions to the Schroedinger equation for this potential are 
gaussians and integrals over gaussians and they may be used to construct 
the Green's functions 

Gi (x, x) = — ^ e"' / e-f' dy f e"^' dz (40) 

Jo Joo 

G2 {x, x) = -e^' / e-y^ dy (41) 



00 



G2 {x, x) - Gi {x, x) ^ ^ e^' / e-y' dy / e""' dz . (42) 

VTl" Joo Joo 
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The sum rule S in this case gives rise to the relation 

poo 

e^' ( 1 - erf{x) f dx = ln2 (43) 

Jo 

where erf (x) is the error function defined by 

erf (x) ^ —= / dy . (44) 

Jo 

This is an interesting integral relation which has arisen from the sum rule 
for the shifted oscillator potential. 



4.3 Quart ic Potential V(x) = 

The low lying eigenvalues for the even states in the quartic potential can be 
obtained from Bender et al (1977) and the eigenvalues for the odd states can 
be extracted from Hioc and Montroll (1975) by appropriate renormalisation 
of the eigenvalues to accomodate the changed strength of the potential. The 
lowest eigenvalues are listed in Table 2. 



n 


^2n 







1.060362 


3.799673 


1 


7.455698 


11.644746 


2 


16.261826 


21.238373 


3 


26.528472 


32.098598 


4 


37.923001 


43.981158 



Table 2: The low lying eigenvalues for the quartic potential 

The eigenvalues for the higher lying states may be calculated using the 
WKB formula in eq.(22) for A'" = 4 to give 

4 

\2n ~ a + (45) 
A2n+i ~ q; ( n + ^ ) (46) 
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where 



a 



6\/7r 



m. 

r(l) 



(47) 



To verify the sum rules the first k terms may be summed exactly using 
the numerical estimates in Table 2 and the remainder may be evaluated by 
approximating by integrals over n and using the WKB estimates which give 



poo 

Jk+h 

noo 



a 

3 
a 



k + 



k + 



Setting k 



4 we get 

51 ~ 0.45003 + 0.31349 

52 ~ 1.22266 + 0.30413 
S ~ 0.77263 - 0.00936 



0.76352 
1.52679 
0.76327 



to be compared with the values from eqs.(29)-(31) calculated using (3 
corresponding to A'" = 4 which give 



= 0.76330 



(48) 
(49) 

(50) 
(51) 
(52) 

_ 1 

6 

(53) 



r(l)r(i) 

in good agreement with the numerical results. The agreement can be im- 
proved by calculating more of the eigenvalues exactly. 



4.4 Particle In a Box N ^ 00 

To check the sum rules for the case A?" — > 00 we proceed as follows. The 
sum rules scale when the strength of the power law potential is changed by 

a factor 7. The scaling factor for the energy eigenvalue is easily found to be 

2 

7^+2 . Therefore if we consider the potential 



Lt 



iV— >oo 



2\x\ 



N 







= 1 



if X < 

if X ^ 
00 if X > 



IT 

2 



(54) 
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which corresponds to a particle in a box with infinite potential at the walls at 
= ^. The eigenvalues for this potential are the squares of integers. Hence 
the sums over the inverses of the eigenvalues can all be found in closed form 
and are given by 

oo 2 

= i:(2„ + 2)- = I- 

n=0 

oo 2 

^2 = E(2^+l)"' = y 

n=0 

S ^ S2 - Si = ^. (55) 
When the scahng factor is taken into account the sum rule in eq.(29) becomes 

^ - [-YJ r (4/3) r (1-/3) ■ 

If the limiting behaviour of the gamma function Ltz^o ~ is used then 
eqs.(56),(30) and (31) give 

S ^ 2Si ^ ^52 = ^ (57) 
in agreement with the results in eq.(55). 



5 General Symmetric Confining Potentials 

In this section we consider general symmetric confining potentials and exam- 
ine the structure of the sum rules. Let and ^2 be two linearly independent 
solutions of the zero energy Schroedinger equation 

(f 

— *(a;) = V(x) ^(x) (58) 

such that 

Lt^-^o 6 = a , Lt^-^o ^2 ^ b X . (59) 

Then a solution $2 which satisfies Lij._»oo^2 — may be constructed as the 
linear combination 

$2 = 6 + c6 (60) 
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and the solution $i may be chosen as either or ^2 to construct the Green's 
functions using eqs.(8)-(10). Then 

Gi{x,y) = + 1^2(x<) $2(x>) (61) 

G2{x,y) = - ^6(a^<) ^2{x>) (62) 

G2{x,y) - G,{x,y) = - ^ *2 (x<) $2 (x>) (63) 

are the Green's functions corresponding to the odd and even eigenstates and 
the difference between them. These Green's functions may be used in eq.(ll) 
to estabhsh the sum rules 

Si = Y.\ = ^ / i2{x)^2{x)dx (64) 

^ \2n+i ah Jo 

^2 = ^— = -— / Clix) ^2{x) dx (65) 



n=0 
00 



J-^ / \n -| poo 

S = i:i^= -- / *iWdx. (66) 



Ar. abc 

n=0 ^ 



By noting that the boundary conditions satisfied by the functions specified 
in eqs.(57) and (58) show that 



d 

dx 

d 



Lt^^o ^ (6 $2) = ahc 
Lt^^o ^ (6 ^2) = ab (67) 



all three sum rules may be written in the compact form 

'dy\ , 1 



dx^ J '"=° A 



/U=o (68) 
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with three different choices of A and / (x) as given by 



A = +S2 , fix) = / (x) $2 (x) dx 



A 



Si, fix) ^ / 6 (X) $2 (X) rfx 



5 




(69) 



If we restore H and the mass parameter of the system then eq.(68) has the ap- 
pearance of a Schrocdingcr Uke equation for the incomplete overlap integrals 
/ (x) evaluated at ,t = with the sums over the inverses of the eigenvalues 
playing the role of the inverse of an energy-like parameter. We have shown 
that the sum over the difference between the inverses of the even and odd 
eigenvalues in a symmetric confining potential can be calculated from the 
integral of the zero energy solution of the Schroedinger equation which is a 
decaying solution in the asymptotic region. 

In this paper it has been shown that for symmetric confining potentials 
it is possible to establish sum rules involving the eigenvalues of odd and even 
states separately. It has been shown that even when the sum rules for even 
and odd states do not converge the difference between them could converge. 
We have shown that for power law potentials simple expressions for the sum 
rules may be derived which show an interesting structure. Eventhough the 
eigenvalues lor N ^ 2 have to be found numerically the sum rules show that 
the inverses of the eigenvalues add up to some simple ratios involving the well 
known gamma functions. This is an interesting result from a mathematical 
point of view. Higher order sum rules such as those emerging from the 
application of eq.(12) indicate even more intricate relationships between the 
solutions of the Schroedinger equation. 
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